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Computer graphics method and device. 



PCT/IB98/D0591 



The invention relates to a computer graphics rendering method according to the 
preamble of Claim 1 and to a device implementing the method. 

Such a method is known from an article by Alyn Rockwood, Kurt Heaton and 

5 Tom Davis, titled "Real-Time Rendering of Trimmed Surfaces" and published in Computer 
Graphics, Volume 23, Number 3, July 1989 pages 107-116, which will be referred to as "the 
Rockwood article" . The Rockwood article describes a way of translating an algebraic model 
of a curved siu^ce patch into a mesh of flat triangles approximating the siurface patch. This 
makes it possible to combine the advantages of the compactness of the algebraic model with 

10 the rendering efficiency of known devices for rendering a mesh of flat triangles. 

The algebraic models used in computer graphics use a parameter space, with 
points having parameter coordinates (u,v). An area from parameter space is mapped to the 
curved surfrice patch, with points having three dimensional coordinates (x,y,z). The 
Rockwood article uses a rectangular area. In the interior of this rectangular area in parameter 

15 space a grid of evenly spaced points is defined. The coordinates of the points on the surface 
patch to which this grid of points maps are computed. In the interior the mesh is constructed 
from rectangles that have these selected points on the surface patch as theu- comer points. 

Points on the boundary edges of die rectangular area are also used to conq)ute 
coordinates of comer points of triangles in the mesh, but the selection of these points on the 

20 boundary edges is treated separately from selection of the grid of points in die interior. The 
reason for this is that a graphics model is usually made up of a patchwork of abutting surface 
patches. In this case, a boundary curve of the surface patch usually coincides with the 
boundary curve of an adjacent surfece patch. To prevent that cracks (or holes) appear 
between the two meshes used to approximate two adjacent surface patches respectively, it 

25 must be ensured that on the shared boundary curve between the adjacent surface patches the 
same points are used for generating both meshes. This is realized by selecting the spacmg of 
the points on each boundary edge of the rectangular area in dependence on the properties of 
the corresponding boundary curve only, so that these points are the same, no matter what 
shape is defined for the rest of the siurface patch. 
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As a result, the spacing of the points on the boundary edges of the rectangular 
area may differ from the spacing of the grid of points facing these boundary edges. This 
means that to form flat triangles a special algorithm is required to connect points of the 
surface patch corresponding to the boundary edge and points of the surface patch 
5 corresponding to the grid. 

Amongst others, it is an object of the invention to provide for a method 
according to the preamble in which connection of points of the surface patch corresponding 
to the boxmdary edge and points of the surface patch corresponding to the grid is simplified 

10 for at least one of the boimdary edges without undue loss of accuracy of the approximation 
of the surface patch by a mesh of flat triangles. 

The method according to the invention is characterized by the characterizing 
part of Claim 1, As a result no special connection measures are needed at least between the 
segmentation points on the selected boundary edge and the grid. It has been found that the 

IS extension of the segmentation step size of a boundary edge to define the grid size of the 
interior generally provides a sufficiently accurate approximation for rendering purpose, 
provided one uses the boundary edge that contains the most segmentation points, or at least 
no fewer segmentation points than any other boundary edge. The approximation is 
particularly good for curved surfru:e patches that have a minimum number of boundary 

20 edges, preferably three boundary edges. 

An embodiment of the method according to die invention is characterized by the 
characterizing part of claim 2. By proper ad^tive selection of the boundary edges between 
whidi the straight lines extend, it is thus ensured that for two of tiie three boundary edges no 
special connection measures are needed. 

25 A iurdier embodiment of the method according to the invention is characterized 

by the characterizing part of claim 3. In a forward difference computation new coordinates 
are computed by adding an offset to preceding coordinates. The offset itself is also be 
computed by adding a further offset to a preceding offtet. For a polynomial of second degree 
in the parameter space coordinates the further offset is constant along the straight line; for 

30 higher order polynomials the further offset is confuted by adding a yet further offset to the 
further offset used for a preceding point along the straight line and so on. The forward 
difference conqputation eliminates the need to perform multiplications to compute the 
coordinates of different points, making it possible to use simple dedicated hardware which is 
able to perform this computation quickly. 
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A further embodiment of the method according to the invention is characterized 
by the characterizing part of claim 4. When the points are evenly spaced along the straight 
line from one boundary edge to the other, computation of the initial values by means of 
forward difference is only possible if the number of pomt on a particular straight line is 
S either the same for eadi particular straight line or increases by a constant number from one 
straight line to the other away from the point where the boundary edges intersect. 

Another embodiment of the method according to the invention is characterized 
in that by the characterizing part of claim 5. At has been found that the curvature of the line 
that corresponds to a boundary edge in the rendered image is the most important criterion for 
10 selecting the required number of segmentation points on that boundary edge. 

Another embodiment of the method according to the invention is characterized 
by the characterizing part of claim 7. By selecting between different methods given the 
number of segmentation points used for the cim^ed surface patch under particular viewing 
circumstances, one may for example minimize the number of generated flat triangles. 

15 

These and other advantageous aspects of the invention will be non-limitatively 
described using the attadied figures, of which 

Figure 1 shows an apparatus for graphics rendering of a modeled object, 
Figure 2 shows the geometry an object described by a number of curved surface 

20 patdies. 

Figure 3 shows surface patches approximated by flat triangles, 
Figure 4 shows a flow chart for generating the set of flat triangles. 
Figure 4a illustrates a m^od of computing a number of segments. 
Figure Sa-e illustrate the generation of a set of flat triangles, 
25 Figure 6a-e illustrate the generation of a set of flat triangles. 

Figure 7a-c illustrate the generation of flat triangles in an interfrure band. 

Figure 1 shows an ^paratus for computer graphics rendering of a modeled 
object. The apparatus contains a memory 10 for parameters that describe the geometry of an 
30 object. The memory 10 is coupled to a transformation unit 12 which transforms the 

parameters. The transformation unit 12 is coupled to a flat triangle generating unit 14 which 
uses the parameters of the model to generate a description of a set of flat triangles which a 
approximate the surface of the object. The triangle generating unit 14 is coupled to a triangle 
rendering unit 16, which computes the pixel values of an image of the flat triangles. Triangle 



wo 98/57301 PCT/IB98AI0591 

4 

rendering units 16 are known per se and are capable of rendering sets of flat triangles very 
efficiently. The triangle rendering unit 16 is coupled to a display unit 18, such as a CRT or 
LCD unit, which displays an image containing a rendering of the flat triangles. At least part 
of the apparatus may be inq)lemented as a suitably progranuned general purpose computer. 
S Parts of the apparatus, such as the triangle rendering unit 16 may also be implemented as fast 
dedicated hardware. 

Figure 2 shows the geometry an object described by a number of curved surface 
patches 20, 22, 24. Curved surface patches 20, 22, 24 can be described for example as 
Bezier triangles. To describe a curved surface patch 20, 22, 24 as a Bezier triangle, the 

10 memory 10 stores the three dimensional coordinates of a nimiber of control points. The 
control points define a mapping from two-dimensional barycentric siu^ce coordinates 
(u,v,w) in a triangular domain (l>u>0,l>v>0,l>w>0, u-hv+w=l) to points on the 
curved surface patch 20, 22, 24. When the model contains neighbouring Bezier triangles, 
these Bezier triangles are defined so that a boundary edge (i.e. either the boundary edge 

15 defined by u=0, or the boundary edge defined by v=0 or the boundary edge defined by 
w-O) of one Bezier triangle coincides with a boundary edge of the neighbouring Bezier 
triangle. For Bezier triangles this is a matter of sharing some of the control points between 
neighbouring Bezier triangles. Thus a continuous description of the surface of the object is 
realized. 

20 When an object described by Bezier triangles is to be rendered, the 

transfonnation unit 12 receives the coordinates of the control pomts and transforms these 
coordinates, for exan4>le to determine the coordinates of the control points in a coordinate 
frame defined by a user selectable viewpoint and view-direction. The transformation umt 12 
supplies the transformed coordinates of die control points to the flat triangle generatmg unit 

25 14. The flat triangle generating unit 14 uses the transformed parameters of the model to 

generate a description of a set of flat triangles which a approximate the surface of the object. 

Figure 3 shows surface patches 30, 32 approximated by flat triangles 34a-i, 36. 
The number of triangles used to approximate the curved surface patches 30, 32 is a 
compromise: the ^proximation will become better as more flat triangles are used to 

30 represent one curved surface patch, and on the other hand the time needed by the triangle 
rendering unit 16 to render the triangles becomes less as fewer flat triangles are used. To 
attain an optimal compromise, the flat triangle generating unit 16 selects the number of 
triangles in dq[>endence of properties of the curved surface patch and circumstances under 
which the curved surface patches 30, 32 are to be rendered: for example in dependence on 
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the curvamre of the surface patch 30, 32, the number of pixels taken up by the curved 
surface patch 30, 32 in the image on the display unit 18, the angle range between the line of 
view to the curved surface patch 30, 32 from the viewpoint, the required image quality etc. 

When the selection of the number of flat triangles for different surface patches 
5 is different for neighbouring surface patches 30, 32, a crack 38 may appear between 
neighbouring surface patches, i.e. in the approximation the surface of the modeled object 
may have areas for which no visual properties are computed. This would lead to undesirable 
artifacts in the displayed image. This is prevented by starting the process of generating flat 
triangles at the boundary edges of the curved surface patch which is to be approximated by 
10 these flat triangles. Each boundary edge is approximated by a set of linked straight line 

segments, which is selected as a function of the properties of the boundary edge, the viewing 
circumstances etc., but independent of properties of the surface patch as far as they don't 
affect the boundary edge. The same selection criteria are used for a shared boimdary edge of 
two adjacent curved surface patches, and therefore the same set of linked straight line 
IS segments will be chosen for the shared boundary edge of those adjacent surface patches. 

Subsequently the mterior of the curved surface patch is ^roximated by flat 
triangles. In this approximation it is ensured that the straight line segments are edges of the 
flat triangles at the boundary of the set of flat triangles used to approximate the interior. 

For a boundary edge of a Bezier triangle (e.g. the boundary edge with u=^0) the 
20 linked straight line segments may be obtained for example by selecting equidistant two- 
dnnensional surface coordmate sets, e.g. (u.v,w)=(0,m/n,l-m/n) where n is a positive 
mteger and m=0,L..n. The Bezier triangle description maps these coordinate sets to three 
dnnensional points along die boundary edge. In ibis case the straight lines connectmg pomts 
m three dimensional space with successive values of "m** are used in the set of straight line 
25 segments. Thus in this case the end*points of the straight line segments are points on the 
boundary edge of the curved surface, which ensures that the straight line segments 
approximate the boundary edge. 

Figure 4 shows a flow chart for generating the set of flat triangles. In a first 
step 40 it is computed for each boundary edge into how many segments the boundary edge 
30 should be split. In a second step 41 a set of flat triangles is computed that contains flat 
triangles which have these segments into which each of the boundary edges are split as 
edges. This is repeated for a number of surface patches. 

The nuniber of straight line segments in that into which a boundary edge is split 
(**n" in the example given above), is preferably selected according to the properties of the 
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edge and the viewing circumstances. In a preferred embodiment the number is selected by 
computing a number indicative of the maximum distance in the rendered image between any 
straight Ime segment used to approximate the boundary edge and the actual (curved) 
boundary edge if that would be shown. This ma ^iniiim distance decreases as the number of 

5 segments increases and the number of segments is selected so that the maximum distance is 
less than a predetermined distance, for example 2 or 3 pixel sizes. 

Figure 4a illustrates a method of computing the number of segments. A straight 
line segment 42 approximating the boundary edge 46 is extended once with its own length 
resulting in an endpoint 43. The location to which this endpoint maps in the image is 

10 computed: if the mzpped segment extends between pi and p2 the endpoint is pe=p24-(p2- 
pl). Also the location in the image where the endpoint 45 pn of the next straight line 
segment 44 maps is computed. The nimiber of segments is selected so that the vector pn-pe 
connecting the location of the two endpoints 43, 45 is less than a predetermined number of 
pixels (e.g. 2 or 3 pixels) long. In case of surface patches described as quadratic Bezier 

15 triangles, the length of the vector pn-pe is ^proximately the same for any segment and 
inversely proportional to the number of segments, which makes it possible to compute the 
desired number of segments from one determination of the distance between the endpoints 
43, 45 for a test value of the number of segments. 

Given the sets of straight line segments used to approximate the boundary edges 

20 of the curved surface patch, a set of flat triangles is generated that approximates the interior 
of the curved surface patch. For this purpose several methods will be described. 

In case each of the boundary edges is approximated by fewer than 4 straight 
line segments a point in the interior of the surface patch is selected, for example the 
barycentre (average of the three comer points), and a set of triangles is generated that have 

25 die selected pomt at thek top and respective straight line segments along the boundary edges 
at their base. In case any one of the boundary edges is approximated by 4 or more straight 
line segments, more comer points for flat triangles are generated for the interior of the 
curved surface patch. 

Figures 5a-e illustrate an embodiment of the generation of a set of flat triangles 

30 in that case. In these and the following figures use will be made of a representation in two- 
. dimensional barycentric surface coordinate space (u,v,w) u+v+w=l. The description of the 
curved surface wiU provide a moping from these two-dimensional barycentric surface 
coordinates (u,v,w) to three dimensional coordinates (x,y,z), as is the case for surface 
patches described as triangular Bezier patches. In actual practice often homogeneous four 
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dimensional coordinates (x,y,z,w) (w=l initially) equivalent to three dimensional coordinates 
are also used to represent a three dimensional point. This is convenient for clipping 
purposes, but does not affect the present invention. Nodes and lines connecting these points 
will be shown in the two-dimensional surface coordinate space. To indicate flat triangles the 
S foUowing convention will be used: when three such nodes are shown connected by lines in a 
triangle in tvro dimensional surface coordinates, this represents a flat triangle having as 
comer points the points to which the nodes map according to the description of the curved 
surface patch only when no further lines are shown crossing the triangle in two dimensional 
surface coordinate space. 
10 Figure 5a shows a surface coordinate space representation of a triangular 

surface patch with boundary edges 50, 52, 54. On each boundary edge 50, 52. 54 nodes 
500a-e, 520a-c, 540a-c have been indicated by dots. Using the mapping from two- 
dimensional barycentric surface coordinates (u,v,w) to three dimensional coordinates (x,y,z) 
these nodes correspond to points in three dimensional space. These points link the straight 
IS line segments used to sqiproxiniate the boundary edges. (In surface coordinate (u,v,w) space 
the boundary edges are represented as straight Imes, as in the case of Bezier triangles, but 
mapped to three dimensional ^ace these boundary edges may become curved). 

In a first step, it is determined which of the boundary edges SO, S2, 54 is 
approximated by the greatest number of segments, or at least by no fewer segments than any 
20 of die odier boundary edges. This boundary edge will be called die first selected edge. In a 
second step it is determined which of die remainmg edges is approximated by the least 
number of segments, or at least by no greater number than the other remaining edge. This 
boundary edge will be called die second selected edge. In figure Sa the surface patdi is 
shown so that the first selected edge SO is shown on the right and the second selected edge 
2S S2 at die bottom. 

The selected edges are used to define nodes in the interior of the triangle. 
Figures 5b-d illustrates the nature of this definition step by step. 

In figure 5b an interface line 56 has been drawn connecting a top node at the 
intersection of the left and right boundary edge to a first node 520a on the left inside the 
30 bottom boundary edge 52. Furthermore horizontal lines have been drawn starting from each 
of the points SOOa-e on the right boundary edge 50 (i.e. from the nodes on the boundary 
edge with the highest number of nodes) to the interface line 56. 

In figure Sc auxiliary lines have been added between the top node and the other 
nodes S20a-c on the bottom boundary edge. Intersections between the auxiliary lines and the 
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horizontal lines are shown. Each of the horizontal lines corresponds to a set of straight line 
segments and the intersections on a horizontal line represent points in surface coordinates that 
are mapped to nodes which link the straight line segments in the set corresponding to that 
horizontal line. Each pair of successive horizontal lines surrounds a band on the surface 

5 patch. Such a band is approximated by flat triangles that have the mtersections as comer 
points. Each intersection on the horizontal lines surrounding the band is used as a comer 
point in at least two triangles, in combination with an intersection on the other horizontal line 
surrounding the band and a next intersection, either on the same horizontal line or the other 
horizontal line. Figure 5d shows an example of how the intersections may be connected to 

10 obtain triangles (note that the connections run along the auxiliary lines, which is a preferred 
embodiment, but that other connections could be used instead of the auxiliary lines also). 

Finally, in figure 5e, triangles are added for the area of the surface patch 
between the boundary edge at the left and the interface line, called the interface band. The 
triangles that approximate the interface band have the nodes 540a-c on the left boundary 54 

15 and the intersections of the horizontal line and the interface line 56 as comer points. 

Preferably, die nodes on the left boundary 500a-e and the horizontal lines are 
spaced equidistandy in surface coordinate (u,v) space. Similarly, the nodes on the bottom 
boundary 520a-c are preferably spaced equidistandy, and therefore the intersections on each 
horizontal line are also equidistant. 

20 Overall, the generation of triangles iUustrated m figures 5a-e involves the 

following steps 

1) compute die re^)ective numbers of segments by which the boundary edges 
should be £q)proximated 

2) selea a first selected boundary edge which is approximated by no fewer 
2S segments than any other edge 

3) select a second selected boundary edge which is approximated by no more 
segments than any other edge 

4) for each node on the first selected boimdary edge: generate the same number 
of evenly spaced nodes along a line (a "horizontal" line in the figure) extending 

30 from that node in parameter space, the line extending in parallel to the second 

selected boundary edge in parameter space, the spacing between the nodes in 
parameter space being proportional to the line number as counted from the node 
opposite die second selected boundary edge 

5) conqiute the coordinates of the points to which the nodes m^ 
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6) form a mesh of triangles connecting the points 

When it is determined that the first selected edge SO (which is approximated by 
no fewer segments than any of die other edges) is approximated by just as many segments as 
a further boundary edge (the left edge), the interface band is not necessary and the horizontal 
5 Imes wiU run from the nodes on the right edge to the nodes on the left edge. In this case, 
each horizontal line contains just as many nodes as the second selected edge (shown as the 
bottom edge 52). 

Preferably, the surface coordinates (u,v,w) map to three dimensional 
coordinates (x,y,z) according to a polynomial fimction, e.g. a quadratic function: 

10 

X = X200 B2oo(u,v,w) + X02O Bo2o(u*v,w) + X002 Boo2(u,v,w) 
+ ^110 Biio(u.v,w) + xioi Bioi(u,v,w) + xoii Boii(u,v.w) 

y = y200 B2oo(u-v,w) + yo2o Bo2o(u,v,w) + yoo2 Boo2(u>v,w) 
+ yiioBiio(u,v,w) + yioi Bioi(u,v,w) + yo^ Bon(u.v,w) 
15 z = Z200 B2oo(u,v,w) + Z020 Bo2o(u,v,w) + Boo2(u,v,w) 

+ zjio Biio(u,v,w) + Zioi Bioi(u.v,w) + zoii Bon(u,v,w) 

Where }^jJ^, y-^^ and z^jj^ (i,j,k=0, 1 or 2) are the coordinates of the control 
points ijk (ij,k=0. 1 or 2); u,v.w with u+v+w=l are barycentric coordinates; Bipj(u,v,w) 
20 = u**i v**j w**k 2/(i!j!k!) (** mdicating exponentiation) are second degree Bernstein 
polynomials. 

The use of a polynomial fimction makes it possible to conqpute the three 
dimensional coordinates for the nodes on each horizontal Ime iteratively, one node after the 
odier on eadi horizontal line, by additions but without requiring multiplication. For example, 
25 using a second degree polynomial, the x-coordinate X| of successive nodes labeUed ''i" 
(i=0,l,2..) on a horizontal line is coiiq)uted using increments "dj" and •e" according to 

di+i = di + e 

= Xi + di+i 

30 

Analogous relations hold for the y and z coordinates y^ z^. To speed up 
processing die computation of the coordinates can be easily implemented in an integrated 
circuit, usmg a dedicated signal processing cucuit. If a polynomial of higher degree is used, 
d wiU be computed by adding a second order increment instead of e, the second order 
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incTCTnent being computed by forward difference like d itself. For realistic graphics rendering 
normals to the curved surface patdi and texture coordinates are needed in combination with 
the x,y and z coordinates for each node. When the parameter space dq)endence of these 
normals and texture coordinates is s^roximated by a polynomial function these parameters 

5 may also be computed using forward differences. Thus the computation of these parameters 
can also be easily inq)lCTiented in an integrated circuit, using similar or identical signal 
processing circuits as for the coordinate conq)utation. 

The initial values Xq do e are readily derived from the polynomial used and the 
parameters of the line along which the nodes lie. For example if the line "m** has constant v- 

10 coordinate v^^ and the u step-size from node to node along the line m is h^, then 

Xo= ^0 Xo02 + Xon Vni(l-V^ 

do= 5^200 ^ + Xoo2(hin^ -2h^(l-Vni)) + ^110 ^n^m + ^lOl^K^ (l-^m)) ' 

xoii K 

15 e= 2h^^(x2oo+%)2-*iOi) 

Of course different initial values occur if a different ^e of polynomial is used. 
The initial values for successive horizontal lines in turn can be computed by forward 
differences derivable from the e3q[)ressions given above, provided that the expressions for the 

20 mitial values Xq^Aq and e are polynomial functions of v^. This is the case if h^^^ and h^^jo. 
are polynomial functions of v^^, for example if h^, is proportional to v^q as is the case in the 
example of figures 5a-e when the nodes and the horizontal lines are evenly spaced and the 
extreme nodes on each horizontal Ime coincide with the boundary edges on the respective 
boundary edges on both sides of the horizontal lines. This forward difference computation 

25 also lends itself to impl»ientation in an integrated circuit using adder circuits. Thus only the 
initial values for the top of the triangle need to be computed and fed to the forward 
difference circuits which can then compute the coordinates of all points for die curved 
surface patch by simple additions. 

h^^ and h^^Voj are also polynomial functions of v^ if h^, has the same constant 

30 value for each Vj^. Figure 6a-e illustrate another embodiment of the generation of a set of flat 
triangles in which the step-size is constant. Figure 6a is similar to figure 5a and shows a 
surface patch in two-dimensional surfoce coordinate space, with boundary edges 60, 62, 62 
mcluding nodes 600a-e, 620a-c, 640a-c. 

In figure 6b the boundary edge 60 with the maximum number ''n" of segments 
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has been selected and two interface lines 66a-b are shown running in parallel (in surface 
coordinate space) with the other boundary edges 62, 64 of the surface patch. The distance 
between an interface edge 66a-b and the boundary edge 62, 64 to which it runs in parallel is 
such that the interface edge would cut I/n^ part from the other boundary edges, if it were 

5 extended to the other edges. Here "n* is the highest of the number of segments used to 
approximate any of the boundary edges. 

In figure 6c horizontal lines are shown in the iimer area enclosed by the 
mterface lines 66a-b and one of the boundary edges 60. The horizontal lines extend from the 
nodes 600a-e on the right boundary edge 60 which contains the highest number of segments 

10 (or at least a not-lower number of segments than any other boundary edge). In 
figure 6d equidistant nodes have been selected on each horizontal line. On each horizontal 
line the leftmost node is on the intersection with the left interface line 66a and the rightmost 
node is on the right boundary edge 60. The number of nodes per horizontal line increases by 
one for successive horizontal lines from top to bottom. The nodes are mapped to points on 

15 the surfece patch in three dimensional space and the mapped nodes constitute the comer 
points of the flat triangles. Zig-zag connections between successive nodes on pairs of 
adjacent horizontal lines are shown, indicating the coimections between the comerpoints of 
the flat triangles. 

In figure 6e a common additional comerpoint 68 is added to both interface lines 
20 66a,b in the middle b^een fte intersection of the two interface lines 66a,b and the 
intersection of their corresponding boundary edges 60, 64. Interface bands between the 
interface lines 66a-b and the boimdary edges 62, 64 have been filled in with triangles that 
have at least one comerpoint on both a boundary edge 60, 64 and an interface line 
66a,b/additional comerpoint 68, the comerpoints being selected from the nodes between the 
25 straight line segments used to approxhnate the boundary edges 60, 62, 64 and the nodes on 
the interfece lines 66a,b. The additional point 68 is added to prevent visual artifacts due to 
the relatively large gap near the intersection of the interface lines, when there are two 
interface lines. 

Interface bands are unnecessary for further boundary edges that are 
30 approximated by the same number of segments as the first selected boundary edge (the 
boundary edge approximated by the greatest number of segments, or at least by no fewer 
segments Aan any other boundary edge). In that case, the extreme nodes on each line 
connecting the further boundary edge and the first selected boundary edge coincide with these 
the first selected boundary edge and the further boundary edge respectively. 
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The method of figure 6a-e differs from that of figure Sa-e only in that the 
spacing between the nodes is the same on each "horizontal" line. Furthermore, any one of 
the boundary edges may be selected as second selected boundary edge (not only the one 
containing the fewest nodes), but if any one of the boundary edges contains as many nodes as 

5 the first selected boundary edge and the remaining boundary edge contains fewer nodes, then 
preferably the remaining edge is used as second selected boundary edge, so that the 
'^horizontal" lines run in parallel to that second selected boundary edge between the two 
other edges that contain mutually the same number of nodes. 

When the nodes are evenly spaced in parameter space, the three dimensional 

10 coordinates (x,y,z) to which the nodes map can again be computed by forward differences, 
using the same techniques as described hereinbefore. The initial values of these forward 
differences for successive horizontal lines in turn can also be computed by forward 
differences because the step-size between the nodes on the horizontal lines in (u,v,w) space is 
constant making h^^^ constant and bj^^Vj^ a linear polynomial in v^^. 

15 A respective interface band needs to be used near each of the boundary edges 

that is ^proximated by fewer segments than the boundary edge that is approxunated by the 
greatest number of segments. It should be noted that diis is in contrast to the method of 
selecting nodes illustrated by figure Sa-e, where no interface band is needed for the lower 
boundary edge 52, even when this lower boundary edge is i^roximated by fewer segments 

20 than the right boundary edge 50. 

Figures 7a-c show the results of several generally applicable methods of 
generating flat triangles in an int^face band. In each example the interface band is shown 
enclosed b^een two lines 70a*c, 72a-c. Each line 70a-c, 72a-c represents a set of linked 
straight line segments in two-dimensional surfoce coordinates. Nodes 700a-c and 720a-c are 

25 shown at the points where the straight Ime segments are linked. The methods shown m 
figures 7a-c make use of the fact that it is easy to generate flat triangles between two lines 
that contain the same number of nodes (plus or minus one) by connecting successive nodes 
from one line to successive lines of the other line and dividing the quadrangles formed in this 
way along diagonals. In the methods the lines 70a-b, 72a-b are divided into parts containing 

30 nodes that can be so connected and into remaining nodes on a line that are connected to a 
conunon node. 

The method shown in figure 7a takes a central part from a first line 70a 
contaming the most nodes, the central part contaming as many nodes as a second line 72a 
containing the least nodes. The nodes of the central part and the of the second line 72a are 
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connected successively. The nodes that remain on the first line 70a to the left and right of the 
central part are connected to the nodes at the endpoints of the second line 72a. 

The method shown in figure 7b divides the second line 72b, which contains the 
least nodes into a left and right part meeting at a central node. Further left and right parts are 

5 taken from right and left extreme of the first line 70b that contains the most nodes. The 
further left and right part contain as many nodes as the left and right part respectively and 
are successive nodes of these parts are connected to successive nodes of the left and right 
part respectively. The nodes that remain on the first line 70b between the further left and 
right parts are connected to the central node on the second line. This method gives a better 

10 image than the method shown in figure 7a if the nimiber of nodes in the two lines 70b, 72b 
differs considerably. 

The method shown in figure 7c is similar to that of figure 7b except that 
crossing diagonals are drawn between the successive connections in the left and right parts. 
At the crossings an additional node is added. Instead of at the intersection of the diagonals, 

IS the additional node may be located at the average of the fom* comers of the quadrangle in 
which it lies» the additional point being connected to the four comers. This method gives a 
better image in case of highlights. 

In case more than one method of generating triangles is available, for exanq)le 
the two methods illustrated in figures Sa-e or 6a-e, a choice is preferably made for each 

20 curved surface patch as to which of the methods is used to approximate it by triangles. This 
is done for example by predicting the number of triangles that needs to be computed for each 
method and using one of those methods that requires no more triangles than the other method 
to generated the triangles. In the case of the two methods described above, the difference 
between the number of triangles can be expressed in Nl, N2 and N3 the number of segments 

25 used for the various boundaiy edges sorted in non-mcreasing order (Nl > :=N2> =N3) 

diff = Nl^ -7N1 +6 if N3< 3 and N2=N1 

diff = Nl^-5Nl+8 if N3<3 and N2<N1 

difi^ = Nl^ -N1(2N3+1)+4N3 if N3> =3 and N2=N1 

30 diff = Nl^ -Nl(2N3-l)-f4N3+3 if N3> =3 and N2<N1 



If the difference "difT is smaller than zero, the method with the same spacing 
between the nodes on each horizontal line should be used, otherwise the method with spacing 
proportional to the distance between the line and the ''top" of the triangle. This choice can be 
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simplified by an iq>proxiinatioii: use proportional spacing if N3<3 or Nl > =2N3-1. To 
inq)rove the quali^ of the approximation this is preferably overruled if the width of the 
triangle in three dimensional space halfway the triangle is greater than at the base: in that 
case constant spacing is preferably used. 
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1. A computer graphics rendering method comprising the steps of 

- defining a curved surface patch in terms of a description of how a polygonal area in a 
parameter space maps onto the curved surface patch; 

- selecting, independently for each boundary edge of the polygonal area, a set of 
5 segmentation points along that boundary edge; 

- selecting a set of grid points inside the polygonal area; 

- computing coordinates for a set of mapped points on the surface patch to which the sets of 
segmentation points and the set of grid points map; 

- rendering an image of a mesh of flat triangles^ nodes of the mesh having the coordinates of 
10 respective points in the set of mapped points; 

characterized in that the method comprises 

- selecting a boundary edge which has no fewer segmentation points than any of the other 
boundary edges of the polygonal area; 

- placing successive rows of the grid points along successive of straight lines in parameter 
15 space extending from successive segmentation points on the selected boundary edge. 

2. A method according to Claun 1 characterized in that the polygonal area is a 
triangle, the method comprising determining whether two boundary edges of the triangle 
have the same number of segmentation pomts, which is no less than the number of 
segmentation points of the third boundary edge of the triangle, and in that case extending the 

20 straight lines between corresponding segmentation points on the two boundaries. 

3. A method according to Claim 2 characterized in that the description provides a 
polynomial relation between coordinates in parameter space and coordinates of points on the 
curved sxuface patch, the grid points being evenly spaced along the straight lines, the 
coordinates of successive grid points along each straight line being determined by forward 

25 difference computation from the coordinates of respective preceding grid points on that 
straight line. 

4. A method according to Claim 3, characterized in that initial values for use in 
said forward difference computation for successive straight lines are computed by a further 
forward difference conqiutation from the initial values for preceding straight lines, respective 
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spacings of the grid points on respective straight line being selected so that the respective 
spacing is equal to a further spacing between the segmentation points on the boundary edges 
between which the straight line extends and nearest grid points on the straight line, a number 
of grid points on each particular straight line being either the same for each straight line or 
5 increasing by a fixed number for successive straight lines away from a comerpoint of the 
triangle where the two boundary edges intersect. 

5. A method according to Claim 1 characterized in that said function is chosen so 
that a deviation between how a segment of the curve between two segmentation points should 
appear in a rendered image and a straight line in the rendered image between the 

10 segmentation points at the ends of the segment is less than a predetermined number of pixels 
of the rendered image. 

6. A method according to Claim 1, characterized in that the polygonal area is a 
triangle, the method comprising selecting a further boundary edge containing no more 
segmentation points than any of the other boundary edges, the straight lines extending in 

IS parallel to the further boundary edge, each straight line containing the same number of grid 
points. 

7. A method according to Claim 1, wherein at least two methods of selecting the 
grid points are available, the graphics rendering method comprising computing which of the 
at least two methods for selecting the grid points generates the fewest triangles, the method 

20 of selecting the grid points that generates the fewest triangles or at least no more than any 
other of the at least two mediod being selected for generating the triangles. 

8. A computer graphics rendering device comprising 

- a memory for storing a definition of a curved surface patch in terms of a description of 
how a polygonal area in a parameter space maps onto the curved surface patch; 

25 - first selecting means for selecting, independently for each boundary edge of the polygonal 
area, a set of segmentation points along that boundary edge; 

- second selecting means for selecting a set of grid points inside the polygonal area; 

- a computation circuit for computing coordinates for a set of mapped points on the surface 
patch to which the sets of segmentation points and the set of grid points map; 

30 - rendering means for rendering an image of a mesh of flat triangles, nodes of the mesh 
having the coordinates of respective points in the set of mapped points; 
characterized in that the device con^rises 

- a selector for selecting a boundary edge which has no fewer segmentation points than any 
of the other boundary edges of the polygonal area; the second selecting means placing 
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successive rows of the grid points along successive of straight lines in parameter space 
extending from successive segmentation points on the selected boundary edge. 

9. A device according to Claim 8, the polygonal area being a triangle, the selector 
detennining whether two boundary edges of the triangle have the same number of 

5 segmentation points, which is no less than the number of segmentation points of the third 
boundary edge of the triangle, the second selecting means in that case extending the straight 
lines between corresponding segmentation points on the two boundaries. 

10. A device according to Claim 9, the description providing a polynomial relation 
between coordinates in parameter space and coordinates of points on the curved siuface 

10 patch, the grid points being evenly spaced along the straight lines, the device comprising a 
forward difference computation circuit for computing the coordinates of successive grid 
points along each straight line being by forward difference computation from the coordinates 
of respective preceding grid points on that straight Ime. 

11. A device according to Claim 10, wherein initial values for use in said forward 
IS diff(^nce computation for successive straight lines are computed by a further forward 

difference computation from the initial values for preceding straight lines, respective spacings 
of the grid points on respective straight line being selected so that the respective spacii^ is 
equal to a further q[>acing betwem the segmentation points on the boundary edges between 
which the straight line extends and nearest grid points on the straight line, a number of grid 
20 points on each particular straight line being either the same for each straight line or 

uu^easing by a fixed number for successive straight lines away from a comerpomt of the 
triangle where die two boundary edges int^sect. 

12. A device according to Qaim 8, wherein said function is chosen so that a 
deviation between how a segment of die curve between two segmentation points should 

25 q)pear in a rendered image and a straight line in the rendered image between the 

segmentation points at the ends of the segment is less than a predetermined number of pixels 
of the rendered image. 

13. A device according to Claim 8, the polygonal area being a triangle, the device 
selecting a further boundary edge containing no more segmentation points than any of the 

30 other boundary edges, the straight lines extending in parallel to the further boundary edge, 
each straight line containing the same number of grid points. 

14. A device according to Claim 8, capable of at least two methods of selecting the 
grid points, ttie device conq>uting which of the at least two methods for selecting the grid 
points generates the fewest triangles, the method of selecting the grid points that generates 
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the fewest triangles or at least no more than any other of the at least two method being 
selected for generating the triangles. 
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